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Abstract 



^ | The photon density operator function is used to describe the prop- 

i— i ' agation of single-photon pulses through a turbulent atmosphere. The 

effects of statistical properties of photon source and the effects of a 
random phase screen on the variance of photon counting are studied. 
A procedure for reducing the total noise is discussed. The physical 
mechanisms responsible for this reduction are explained. 



1 Introduction 

Fluctuations of the atmospheric refractive index caused by turbulent eddies 
considerably limit the performance of long-distance optical communication 
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systems. An initially coherent laser beam acquires some properties of Gaus- 
sian statistics in the case of a long propagation path or strong turbulence. 
The tendency of the scintillation index (defined as a normalized intensity 
variance) to approach asymptotically the level of unity in this case is a dis- 
tinct manifestation of gaussianity of the field statistics. In the literature this 
effect, known first from experimental studies of Gracheva and Gurvich [I], 
is referred to as the saturation of fluctuations. (See, for example, [2]- [7].) 
In that way, the strong influence of the refractive index fluctuations on the 
radiation field statistics is emphasized. 

It has been shown quite generally [8] -[10] that the normalized inten- 
sity variance approaches unity for any source distribution as the turbulence 
strength tends to infinity. This result is valid for any degree of coherence of 
the source, provided the response time of the recording instrument is short 
compared with the source coherence time (fast detector). In the opposite 
case of a slow detector, the scintillation index can decrease for partially co- 
herent beams. The effect is not very strong for temporal partial coherence 
11], [H]-[I3], but is very pronounced in the case of spatial partial coherence. 
(See papers [Hj-[T6]. which deal with stationary beams.) 

The purpose of the present paper is to study the effects of partial spa- 
tial coherence on the statistics of the detected photons when photons are 
generated as individual pulses of electromagnetic field propagating in the 
earths atmosphere. Moreover, the effects of the initial statistics of photons 
on fluctuations of the detector counts will be elucidated. 

The case of single-photon pulses is of special interest for quantum cryptog- 
raphy because individual photons are carriers of information bits in several 
basic strategies for free-space quantum key distribution. (Practical free-space 
quantum key distribution is described in references [I7]-[IH].) Fluctuations 
of the detector counts for the case of single-photon pulses were studied both 
theoretically and experimentally in [20]. In contrast to [20] where two limit- 
ing cases (plane-wave and spherical-wave approximations) were studied, our 
consideration takes into account the actual beam diameter. Also we will not 
restrict ourselves with the case of small fluctuations of the radiation field as 
in [20] • The case of strong fluctuations will be studied as well. The method 
of photon distribution function in phase space developed in [16] is somewhat 
generalized here to apply to the case of laser pulses. 
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2 The photon distribution function, and pulse 
propagation in the atmosphere 

Similar to |16j . we proceed from the quantum version of Hamiltonian, H, of 
photons in a medium with a fluctuating refractive index, n(r) [n(r) — 1 << 1] 

h = J2 fi^kK bk-Yl ^k^k' K & k+k' • ( i ) 

k k,k' 

Here the terms describing the zero-point electromagnetic energy are omitted; 
the two terms on the right-hand side describe photons in a vacuum and the 
effect of refractive index fluctuations, respectively; and bk are creation 
and annihilation operators of photons with momentum k, Tiuo^ = hck is the 
photon energy; c is the speed of light in a vacuum, and n-^ is the Fourier 
transform of n(r) defined by 

n k = ^ / dV^nir), (2) 

where V = L x L y L z is the normalizing volume. 

Eq. ([T]) is valid in the limit of small wave- vectors k' (k' <C k). This means 
that the scale of spatial inhomogeneity of turbulence is much greater than the 
wavelength of the radiation. Also we assume here that the initial polarization 
of light remains unaffected by the turbulence throughout the distance of 
propagation. The depolarization of light due to atmosphere turbulence is 
very small. (See, for example, [21] and [22] where this effect was studied.) 

The photon distribution function is defined by 

/(r,q,t) = ^E^ kr < +k/2 W- (3) 
v k 

The operator function /(r, q, t) is the photon density in six-dimensional 
(r, q) phase space. We will use it to describe light pulses with characteristic 
sizes, I, much greater than the wavelength of the radiation, A. In this case, it 
is convenient to restrict the sum over k by some k (k < k « q , k > 2n/l, 
where go is the wave vector corresponding to the central frequency c^o of 
radiation, uj = cqo). Then, the kinetic equation for the distribution function 
is governed by 

{3 t + c q ^ + F(r)0 q }/(r,q,t) = O > (4) 



3 



where F(r) = Uod r n(r) an d c q = du q /dq. When deriving Eq. (j3J) we have 
considered the refractive index to be a slowly varying function of the coordi- 
nate, r. (See more details in [16].) In this case the effect of the turbulence 
on / is represented by the random force F(r), that can be seen in Eq. 
The general solution of Eq. (TjJ is 

n ^, t) ^{r-jy^-jy-fi ;t0 }, ( 5) 

where the functions r(t') and q(t') are photon "trajectories" defined by the 
equations of motion 

Sl(t) «[,(*)], 



dt 

0q(t) 



F[r(t)], (6) 



m 

and the corresponding initial conditions. The last equations follow from the 
requirement that the trajectories in Eq. (|3j) r(t') and q(t') pass through the 
point r, q at t' — t [i.e. r(t' = t) = r, q(t' = t) = q)]. The initial value of 
/(r,q,t), is 0(r,q,t o ) 

0(r,q,t o ) = /(r,q,t ) = tt E e" ikr (6j +k/2 Vk/2)| t = to = E e" 4kr 0(k, q, t ). 

v k k 

(7) 

It is convenient to set to equal to the instant just after photon exits 
the source (as shown in Fig. 1). These photons have not been affected by 
the atmospheric turbulence and their statistics is determined by the source 
properties only. The operators, b*, and 6 q , describe amplitudes of the field, 
E atm , outgoing from the source. The correspondence between free-space and 
generated modes can be established using the following reasoning. The field 
in the atmosphere is given by 

E atm (r) = i E i^y^j V ' V qr & q - e~^b+). (8) 

On the other hand, the outgoing field localized in the vicinity of the laser 
aperture can be expressed in terms of the laser mode creation and annihi- 
lation operators, b + and b, and the normalized functions, $*(r) and $(r), 
describing the spatial distribution of the field, as in Ref . [23J : 

E s (r) = i(2Trhtu ) 1/2 [exp(iq z)^(v)b - exp(-iq z)<f>*{r)b + ). (9) 
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Figure 1: The scheme of the communication channel. 



The function, $(r), can be chosen as 




where r\ = x 2 + y 2 ; z is the position of the pulse center at t = t ] and r 
and r z describe the aperture radius and the length of light pulse, respectively. 
Then, from the condition E atm (r,t Q ) = E s (r,t Q ) we obtain 

6 q (t ) = ^7=/ dre-^-^ r $(r). (11) 

It should be noted that the condition E atm = E s is not rigorous. We ignore 
here vacuum fields arising from other modes of the resonator as well as the 
zero-point free-space waves reflected from the output window. These fields, 
providing correct commutation relations for the operators &+, 6 q , do not 
contribute to the detected number of photons. There are no photons in these 
fields. (Direct detection is assumed.) Thus, the corresponding "vacuum" 
terms are omitted in Eq. (II ip as being irrelevant to our problem. 

The effects of the phase screen can be included by introducing the Mul- 
tiplier, e iip ^ T± \ into the integrand of Eq. ffTTj) . where (p(r±) = ar_i_ and a is a 
Gaussian random variable with a covariance ((a x>y ) 2 ) = 2A~ 2 . This results in 
the following relation for the average of the exponent (e iar± ) = e~ r ± Xc . This 
relation will be useful in our further analysis. Also, the fluctuations of the 
refractive index are usually considered as Gaussian variables with a spatial 
Fourier- component of the correlation function (n(r)n(r')) g = ^(g) given by 



(12) 



Eq. f)12p is refereed to as the von Karman spectrum. L and Iq are the outer 
and inner scales sizes of the turbulent eddies, respectively. In atmospheric 
turbulence, Lq can range from 1 to 100 meters, and Zo is of the order of 
several millimeters. C 2 is the index-of- refract ion structure constant. In most 
physically important cases the quantity Lq 2 in the denominator of Eq. (j!2p 
can be omitted. In this case, the von Karman spectrum is reduced to the 
Tatarskii spectrum [2J. 

Eqs. ([5]) and (JjO, [T0irT2l are sufficient to determine the beam intensity 



(J(r,t)) = c^^ q (/(r,q,t)) 



(13) 



at any r and t. Represent Eq © in the form 



/(r,q,t) = E ex Pi - «k{r - c q (t - t ) + - f dt'{t' - t )F ± [r(t')]} 1 x 
k t <?0 J to ) 



^ k |q-£^'F ± [r(0];to|. 



(14) 



This form is more convenient for obtaining the explicit form of Eq. (fl~3|) . 
We have neglected the changes of the longitudinal photon momentum (q z ) 
caused by the turbulence, because of their negligible contribution to the 
initial momentum (sa q ) for almost any reasonable propagation path. 

The photon "trajectories" r(t') in the arguments of can be approxi- 
mated by straight lines, r(t') = r + c q (i' — t). Thus, we disregard the variation 
of photon momentum at distances of the order of the turbulence correlation 
length, L , which is assumed to be much less than the total propagation 
path, c{t — t ). Then, the explicit expression for (/) is given by 



UM)> 



2 chujc) 



ix nR 2 r z 



exp 



ri-2^ 



R 2 



(15) 



where 



z eff = z -zo-c(t- t ), R 2 = ^U + 



2c{t - 1 ) 



+ 



k 3 (t-t ) 3 r ' 



l + 2r 2 A- 2 ' 
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and 

T = 0.558C^ 1/3 . 

z e ff and R are the distance to the pulse center at time t and the beam radius, 
respectively. 

The effects of partial coherence on the beam radius is represented by 
the quantity, r±, which enters the second term in expression for R 2 . This 
term describes the diffraction broadening of the beam in the course of its 
propagation. The third term is due to turbulence. It dominates at large 
distances. It does not depend on the initial partial coherence. Therefore, the 
beam radius is not sensitive to the presence of phase screen when c(t — t ) — > 
oo. In this connection an important question for our analysis arises: Is there 
any effect of the partial coherence on the photon-counts statistics at large 
propagation paths? This issue is considered in the next section. 

Eq. ffT5]) is derived for a single-photon pulse, i.e. for (b + b) = 1. For arbi- 
trary photon number per pulse, N pu i se , the coefficient in front of the exponent 
in Eq. (TTSJ) should be multiplied by N pu i se . The case of a homogeneous beam 
(see, for example, [TB]) can be considered by assuming r z — > oo and by renor- 
malizing the coefficient in front of the exponent. One can easily see that the 
beam radius is increased with the distance (or propagation time (t — to)) 
similar to the case of the stationary intensity. In contrast, the pulse length 
remains unaffected by the turbulence, within the formalism described above. 



3 Fluctuations of photon counting 

As before, we assume that the background radiation noise is negligible. Also, 
the detector area is taken to be small compared with the beam width. This is 
quite reasonable for systems with long transmission distances. The counting 
interval, T p , is much greater than the pulse duration r p {r p ~ r z /c, 
T p » t p ), as shown in Fig. 2. (In experiments [20], the radiation was in 
the form of Ins laser pulses with a repetition rate 1 MHz and the interval 
T p = 10ms.) 

We start with the definition of the mean-square of the photon counts, n, 
for the time interval, T p . Usually, this is represented as a sum of two terms 

(n 2 ) = (n) + (n(n- 1)), (16) 

where (n) = aN, and a describes the light collection and the quantum detec- 
tion efficiency. It is the ratio of the numbers of detected (n) and generated 
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Figure 2: A sequence of single-photon pulses. 



photons, N. The average number (n) is determined by the integral of light 
intensity I(r, t) as 

rt+T/2 

(n)= V dt'(I(r,t')), (17) 

Jt-T/2 

where rj describes the detector efficiency By substituting Eq. ([15]) into Eq. 
(jT7|) . and assuming the detector is at the beam center, rj_ = 0, we have 



(n) = rjN 



71 



(18) 



Thus, we see that a and r\ are related by 



a 



V- 
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The second term in Eq. f|T6l) is determined by [2 

rt+T p /2 



(n(n-l)) = (:{ V J t Tp P /2 dt'l(r,t')} 2 : ), 



(19) 



where the symbol (: ... :) indicates normal ordering of the operators b^(t') 
and bqi{t") with subsequent averaging. In our case these operators enter the 
definition of the intensity I{r,t') and, hence, the right part of Eq. (fT9|) . 

Each integration over t' within the total interval T v is reduced to a sum 
of N independent integrations within much smaller intervals At, [(r z /c) << 
Ati << T p /N,i = 1,2, ..N], which correspond to instants when the pulses 
cross the detector plane. Then Eq. fTl9|) is reduced to 



(n(n - 1)) = r] 2 N{N 



-lrf^YEE Jdt' Jdt"e-^ +k '>: 

\ I q,k q',k' 
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(^Va(06} +W2 (O6q'-v/a(OW(0>. (20) 

where integrations are within any two different intervals Aij. The coefficient 
N(N — 1) indicates the number of such intervals. Coinciding intervals do not 
contribute to Eq. (fl9]) because of the zero value of the intensity-intensity 
correlation for a single-photon pulse. 

Not all terms in the sum over q, k, q', k' contribute significantly to 
{n(n — 1)). The analysis in [IB] shows that at large but finite propagation 
distances the terms with small values of (i) k, k' or (ii) |q' — q + (k + 
k')/2|, |q — q' + (k + k')/2| (almost diagonal terms) are the most important. 
The evolution of such terms can be described in a manner similar to the case 
of the evolution of /. Then, the explicit form of the expression (|2"0"|) is given 
by 

(n(n - 1)) = V 2 N(N - 1)(2tt) 5 (^|^) 2 £ £ e -[Ad+Arf+«-Hfc?)/4]^/2 x 

* q,k q',k' 

5(fc z )5(^)e-( Q2+Q ' 2 )'-?/ 2 -^ +fe 2)^/8 + 

5(g 2 _ q' z )6(k z + ^) e -[(Q+Q') 2 + (kx-kY) 2 /4]r-?/4-[(Q-Q') 2 +(k x +k' ± ) 2 /4]r2/4 
e - 4 {k x [r-c(q)t,+^ ^ dt 1 t 1 F(r q (t 1 ))]+k' x [r-c(q%+£ /„*' *i*iP( V (ti))] }\ ^ 

where Aq z = q z -q ; Q = q ± -/ ' z dt 1 F ± (T^(t 1 )); Q' = qWo* ^iF±(r q *(ti)), 
and the function r q (ti) is the particle trajectory which passes through the 
point (r, q) at the instant t\ = t z = t — t . 

Summation over the variables q z , q' z , k z , k' z can be easily performed. 
After that, the remaining sum coincides exactly with the corresponding sum 
in the right part of Eq. (33) of reference [16J . Hence, we have 

(n(n-l)} = a 2 N(N + a 2 ), (22) 

where a 2 is the scintillation index for the stationary beams calculated in 
[T6] using the formalism of photon distribution function that is similar to 
the presented here. It follows from Eqs. (120]) and (1221) that the normalized 
variance of photon counting can be written in the form 

W = TT + I n- (23) 
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Eq. (I2"3"j) was obtained for the case of long propagation distances. Neverthe- 
less, it has the same form in the opposite limit of short distances (or weak 
turbulence) when perturbation methods like Rytov's approach [2], [16] are 
applicable. 

The first term in the right part of Eq. (123]) is due to the discrete (quan- 
tum) nature of the photon field. When a — > 0, it is reduced to 1/ (n) which 
is known in the literature as the shot-noise limit (or the standard quantum 
limit). The coefficient (1 — a) in the numerator arises from the (n(n — 1)) 
term in Eq. ({TBI) and is evidence of the nonclassicality of the light. The pres- 
ence of the factor (1 — a) in Eq. (123]) is typical for squeezed photon- number 
radiation. In our case, the total photon number, N, is considered to be a 
constant. For photocount statistics, this situation is equivalent to the case 
of photon number state (Fock state) of the light. In the hypothetical case of 
a — 1, the quantum term vanishes. This case corresponds to the physical sit- 
uation in which all generated photons reach the detector aperture and each 
photon produces a photocount. Of course, this situation is impossible for 
long propagation distances because the beam radius becomes much greater 
than the receiver aperture. In addition, the detection quantum efficiency is 
always less than 100%. 

The second term in Eq. (|23|) can be interpreted as being caused by at- 
mospheric turbulence. The scintillation index determines its relative contri- 
bution to the total noise. In the general case, a 2 depends on the propagation 
distance, the radius of the source aperture, turbulence strength, and the cor- 
relation length, A c . Fig. 3 shows that the value of a 2 can be suppressed 
considerably by decreasing the initial coherence length, A c . Significant ef- 
fect is achieved even for a moderate decrease of the dimensionless parameter 
( r i/ r o) 2 - F° r example, we see in Fig. 3 an almost 50% decrease in a 2 when 
(ri/ro) 2 = 1/2. 

When deriving Eqs. (I2"2l [231) . we have assumed that the total number of 
photons N in the interval T p is constant (no fluctuations of N) . In general, N 
is a fluctuating quantity. For example, much research uses heavily attenuated 
laser pulses to approximate a single-photon source [J7] , [J5] , [2U] , [25] . (Usu- 
ally the photon number per pulse is less than 1.) In this case, the random 
variable iV obeys Poisson statistics resulting in the average value (N(N — 1)) 
equal to (N) 2 . Then, Eq. ( 1221 becomes 

(n(n-l)} = a 2 (N) 2 (l + a 2 ). (24) 
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Figure 3: Dependence of scintillation index on the turbulence strength G\. 



Instead of Eq. fl23|) we have 

W = V) + {) 

Comparison of Eqs. (123]) and (1251) shows the increase of total noise in the 
latter case that is due to contribution of the generation-rate fluctuations. 



4 Conclusion 

The central point of the present paper is the possibility of reducing the photo- 
count fluctuations. It follows from the previous Section that a single-photon- 
on-demand source or squeezed photon-number light provide lower noise levels 
than a heavily attenuated laser source. Also, an additional decrease in pho- 
tocount fluctuations can be achieved by means of a random phase screen. 

It can be seen that the phase screen can decrease the count noise for long- 
distance propagation. In the opposite case a similar effect can be achieved 
by improving the noise characteristics of the source. 
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The question arises, "What is the physical nature of the second terms in 
Eqs. (I2"3"j) and (I2"5"j) for the case of single-photon pulses?" These terms de- 
scribe intensity-intensity correlations. At the same time, one can easily see 
that this correlation is absent within any given pulse. Also, different pulses 
are independent events. In principle, they can be generated even by different 
sources. Different photons do not interact one with another. Nevertheless, 
a 2 7^ 0. This paradox can be explained as follows. Different photons move in 
a fixed distribution of the refractive index. (It was tacitely assumed that the 
turbulence does not vary during the integration time T p .) They are affected 
by the same random force, F. Therefore, there is a strong correlation of 
photon's trajectories. The probabilities of detection for any photons prop- 
agating within small time interval T p depend on the same refractive- index 
configuration. Hence, the detection events within this interval T p are corre- 
lated. The main purpose of the phase screen is to destroy this correlation. 
When the characteristic time of phase variation introduced by the phase 
screen is of the order r d << T p and, in addition, A c << r , the scintillation 
index is decreased considerably for long propagation paths. Thus, a rapid 
spatio-temporal variation of the phase is an effective method for decreasing 
the photon-counting noise. 
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